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ABSTRACT 

A static configuration of point charges held together by the gravitational 
attraction is known to be given by the Majumdar-Papapetrou solution in the 
Einstein-Maxwell theory. We consider a generalization of this solution to non- 
Abelian monopoles of the Yang-Mills Higgs system coupled to gravity. The so- 
lution is governed by an analog of the Bogomol'nyi equations that had played a 
central role in the analysis of non-Abelian monopoles. 



Under the combined forces of Coulomb repulsion and Newtonian attraction, one can 
envisage static configurations of any number of point masses at arbitrary locations. This 
happens when the charges, all of the same sign, are equal in magnitude with the correspond- 
ing masses in certain units so that the Newtonian attraction is balanced by the Coulomb 
repulsion. It is remarkable that such a static configuration is allowed in the framework of 
general relativity as a solution of the Einstein-Maxwell equations. This was first shown by 
WeylJIJ assuming axial symmetry and later generalized by Majumdar and Papapetrou[^]. 
Hartle and Hawking have interpreted this solution as representing an assemblage of ex- 
treme Reissner-Nordstrom black holes. Due to the inherent symmetry between the electric 
and magnetic fields in the Maxwell theory, one can replace the electrically charged point 
masses of the Majumdar-Papapetrou solution by magnetic monopoles. In some respects, 
this suggests an analogy with the case of magnetic monopoles in the framework of sponta- 
neously broken non-Abelian gauge theories. When the Higgs field responsible for breaking 
the symmetry belongs to the adjoint representation of the gauge group, there exists a 'no- 
interaction' result [f|] for a system of monopoles governed by the Bogomol'nyi equations]^] in 
the Prasad-Sommerfield limit Q. One is led to a very simple physical interpretation of this 
result as due to the balance of Coulomb-type repulsive forces among the magnetic charges 
and Yukawa-type attractive forces due to Higgs fields. A natural question that arises is what 
happens when such monopoles are coupled to gravity. In this letter, we consider a type 
of coupling to gravity that admits static solutions that may be viewed as a generalization 
of the Majumdar-Papapetrou solution to a non-Abelian theory or as an extension of the 
Bogomol'nyi equations to include gravity. These solutions represent static configurations 
of non-Abelian monopoles held in place by the gravitational forces rather than the Yukawa 
forces. 

The coupling to gravity of the Yang-Mills Higgs system we consider is given by the 
following action J7J: 

s = / A^-^-^Rf- - i ( F *„y + i w -\(4?- «>)* } , ID 

where R is the Ricci scalar, F* v is the field strength associated with the gauge field and 
(f) a is the Higgs field in the adjoint representation. -D M a is the gauge covariant derivative of 
<p a and (p 2 is a short form for <j) a (j) a . When <p 2 = v 2 , the first term gives the usual Einstein 
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Lagrangian. The rest of the terms represent the standard Yang-Mills Higgs system. We seek 
a static solution in this system assuming Aq = 0. We will set our unit of scale to correspond 

tO 4:7fG = 1. 

In the absence of gravity, one can obtain a static solution in closed form in the Prasad- 
Sommerfield limit (A — > 0) by use of Bogomol'nyi equations. In this limit, in terms of the 
'electric' and magnetic fields 



E? = A0 a , Bf = -e ijk F? k , (2) 



one deduces the following inequality for the energy functional: 
£ = \jd?x[{E?f + {Bt) 



= ^Jd 3 x (E? t Bf) (Ef t Bf) ± J d 3 xE?B? 

> ± / d 3 xE?Bf, (3) 



where the last integral can be written as a surface integral using (g) and the Bianchi identity 
for Fij. Clearly, this becomes an equality only if the Bogomol'nyi equations Ef = ±5" are 
satisfied. For the gauge group SU(2), a spherically symmetric solution to these equations 
representing the 't Hooft-Polyakov monopole (or antimonopole if — sign is chosen for ±) has 
been found in closed form[[5], ||. 

We find it remarkable that this system, when coupled to gravity as in ([]]), still admits 
monopole solutions with an ansatz that resembles the Bogomol'nyi equations. This happens 
when the dimensionless parameter 4tcGv 2 becomes unity, i.e. in our units v — 1. First, 
motivated by the Majumdar-Papapetrou solution, we make the following ansatz for the 
metric: 

ds 2 = g^dx^da? = V 2 dt 2 - — dzW, (4) 

where V, independent of t, is a function of x l 's. The components of the Ricci curvature are 
then easily computed: 



R 00 = V 3 d 2 V-V 2 (dV) 2 , R 0l = 

d 2 v jdv) 2 

V V 2 



Rij 



*,-2^, (5) 
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where d refers to differentiation along the spatial directions. This yields for the Ricci scalar 
the following expression: 

R = -2Vd 2 V + 4 (dV) 2 = 2V 3 d 2 (—) . (6) 
In the background of the above metric, we wish to minimize the 'energy functional' 

£ = fd 3 x {^- 2 R<P 2 + \v*F*F* + ^A0 a A0°) , (T) 

obtained from the action ([!]) ignoring the A— term for the moment. The 'electric' and mag- 
netic fields are defined to be 



1 V 

rpa r-> v a r>a 77a 



where X a = V<p a . It is now straightforward to rewrite the above energy functional in terms 
of E and -B||. The result agrees with the previous expression in Eq. @ and hence obeys 
the same inequality. The lower bound is still given by a surface integral since the factors of 
V coming from the above definitions cancel out in the product E^Bf. Again, the equality 
holds only if the equations = ±-B? are satisfied. In the present case, this implies 

F$ = ±^e m D k X a . (9) 

This resembles the well studied Bogomol'nyi equations of flat space theories except for the 
factor 1/V 2 . 

Since the above procedure yields a minimum of the energy functional, the Euler-Lagrange 
equations of motion for A ^ and <p a will be automatically satisfied. It follows from ([!]) that 
these equations are 

2 rria\ „ rabcib 



D t {V'F%) = gf^DjP, (10) 
V 2 Di(DiF) = ^R<t> a + \(<t> 2 -v 2 )<t> a , (11) 

where g (not to be confused with the determinant of the metric) is the gauge coupling 
constant and f abc 's are the structure constants of the group. It is now easily checked, using 
that Eq. fllCf ) is identically satisfied. The ansatz (g) requires, as a consequence of the 
Bianchi identity for i 7 ^-, that X satisfy 

D * (^ Ar ) = °- ( 12 ) 



Using this we find that Eq. (|TT|) is also satisfied (again, momentarily ignoring the A— term) 
when the Ricci scalar is given by our previous expression in Eq. (|6|). Thus the equations of 
motion for the matter fields are all automatically satisfied. 

The gravitational equations are obtained by varying the action ([!]) with respect to the 
metric g^ v . They take the form 



1 2 

R^U ~~ ~^9llvR = ~~~jP^lLVl (13) 



where T„ v is given by 



T, u = V 2 F^ + D,rD u <p a + g, u ^F^F* + -V 2 D^D^ + - (0 2 - v 2 ) \ , (14) 



\V*F«F« + Ul 2 D^D^ + \ 

upto terms involving derivatives of <fi 2 . Remarkably, these are all satisfied when the following 
holds: 

(f) 2 = v 2 = 1, or X 2 = V 2 . (15) 

This can be easily verified using Eqs. (Q), (||) and (|i~2"D. We note that our neglect of the 
A— term is justified since it vanishes under this condition. All that is required of the scalar 
potential is that it and its first derivative with respect to <fi 2 vanish at <ft 2 = v 2 . 

It remains to solve Eq. (Q), or flT2"|), to obtain the final solution. Upto now we have 
not made any assumption about the gauge group. If we have U(l), we can drop the indices 
a, b, ■ ■ ■. Then X 2 = V 2 implies X = ±V and hence from Eq. ([[D we get d 2 (l/V) = 0. The 
resulting metric agrees with the well known Majumdar-Papapetrou solution for the Einstein- 
Maxwell theory 0. In this case, the scalar field is completely frozen at its vacuum value. The 
solution corresponds to a collection of point monopoles held together by the gravitational 
attraction. 

Also, we have not made any assumption about the number of monopoles or the symmetry 
of the solution. There should exist solutions representing static configurations of non-Abelian 
monopoles, but unlike those of flat space theories, such configurations will be held in place 
by the gravitational forces rather than the Yukawa forces. This is because the Higgs field is 
effectively at its vacuum value as can be seen from (|T5|). Balance of forces requires that any 
such solution has, in our units, its total mass equal to its total monopole charge. Note that 
the surviving magnetic field at spatial infinity is 

B t = = ±^- 2 D l X a X a = ±^- 2 d t (X 2 ) = ±i$V, (16) 

4 



where we used Eqs. (JD, (||) and (|T5|). Now, at spatial infinity, V is expected to have the 
behaviour 

V r^oc, (17) 



where r is the radial distance and M is the total mass of the system. Using this in ([16]) 
we obtain the total monopole charge, given by the total magnetic flux emanating from the 
system, to be equal to ±M that agrees in magnitude with the total mass. One can also get 
this result from the minimum of the the energy functional: 

£ = ±J d 3 xE?B? = J d 3 xE?E? = J d 3 x^DiX a DiX a , (18) 

which after partial integration using fljjD reduces to a surface integral to give the total 
monopole charge. 

Next we consider the gauge group SU(2) and specialize to the case of spherical symmetry. 
As in the case of the 't Hooft-Polyakov monopole, we start with the following ansatz for X a 
and A\: 

X a = Tr a V(r), A"= — e aib f b (l-W(r)), (19) 

gr 

where f is the unit vector along the radial direction. This gives, for the electric and magnetic 
fields defined in (§), 

r V V r J 

Bf = -^—- — (rW'+l-W*)^ (20) 

where prime denotes differentiation with respect to r. The Bogomol'nyi equations Ef = ±Bf 
then lead to the following nonlinear system: 

r 4 F" = e" 2F , (21) 

?- = - + (22) 
V r 



where F is defined by W = e~ F jr. Eq. (pl|) is a special case of the so-called Lane-Emden 
equations that play a central role in the internal constitution of stars0. To our knowledge, it 
has not been solved analytically. However, one can show that there exists a class of solutions 
with the asymptotic behaviour 

F —> Cr + D , r -> oo, (23) 

F -> -lnr , r -> 0, (24) 
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where C > and D are two constants. The expected behaviour of V given by (|TTD thus 
follows from (|22|) and ([23|) if C = g and M = Air/g. This is in accordance with our earlier 
result that the mass and the monopole charge coincide. However, the behaviour of V as 
r — > is not that simple. Direct use of the asymptotic behaviour of F gives V = oo. We 
need the next correction to this result to come to any conclusion. It is easily found that the 
correction to Q23| ) is of the form 

Ay/r sin ^-y-hrr + , (25) 

for some constants A and 5. Use of this suggests that V goes to zero as y/r, as r — > 0. 
However, it does so in an oscillatory way. The area of any sphere surrounding the origin, 
Airr 2 /V 2 , vanishes periodically and tends to zero in this limit. The Ricci scalar given by 
(§) goes as — 4V 2 /r 2 and diverges whenever the area vanishes, the geometry being singular 
for those values of r. What is then the nature of the singularity at r = ? In the case 
of the Majumdar-Papapetrou solution for the metric, the area tends to a constant and one 
approaches the event horizon of the black hole in this limit 0]. In our case the origin does 
have an event horizon, but with zero surface area. This suggests that the singularity at the 
origin is close to being naked. There are no other horizons to dress this singularity since 
1/V behaves smoothly, and hence V does not vanish, for finite values of r. 

It seems that the above features are specific to AkGv 2 being unity and to get a better 
understanding, one needs move away from this critical value. We may test the sensitivity of 
our conclusions by modifying the equation for F as follows: 

r^F" = e" 27F , (26) 

where 7 is somehow related to AttGv 2 and when unity yields the critical case discussed earlier. 
It is now easy to note that the asymptotic behaviour of F as r — > 00 is unchanged when 7 
is greater than 1/2, but for r — > one finds 

F^-( 2 -i)lnr-i-ln(2-i), r -> 0. (27) 

This shows that near origin V approaches zero as r, 

V -> -^—r, r^O. (28) 
1-7 

6 



This is what is expected for a black hole with a finite event horizon at the origin. This 
sounds attractive, but we have not been able to derive this case from an action principle. 

The constant D of Eq. fl23|) is not determined by the data at spatial infinity, leading to a 
one parameter family of solutions. It appears to be a measure of the inverse size of the core 
containing non-Abelian structure. As it tends to oo the core collapses to the origin and the 
solution coincides with that of Majumdar and Papapetrou for a single monopole. This case, 
where W is zero and 1/V is 1 + 1/gr, can be verified to follow directly from Eqs. (|20p . It 
lies entirely within the U(l) subgroup of SU(2) defined by the direction of the Higgs field. 

To summarize, we have analyzed a type of gravitational coupling to non-Abelian monopoles 
and obtained a solution that may be viewed as a generalization of the Majumdar-Papapetrou 
solution or as an extension of the Bogomol'nyi equations. The case of conventional coupling 
of gravity to non-Abelian monopoles have been discussed in the literature [0] drawing sup- 
port from numerical treatments. The advantage of the present approach is that it admits 
analytical treatment to a large extent. Perhaps, this is possible only in the critical case 
discussed in this letter and moving away from criticality may require numerical analysis. 
Secondly, the present approach opens up the possibility of analyzing multimonopole solu- 
tions in the presence of gravity and it will be interesting if such solutions exist and the 
known techniques are applicable in analyzing them. It will also be interesting if there exists 
a stationary generalization analogous to the one that exists for the Majumdar-Papapetrou 
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